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Confirmation of the y(4260) Resonance Production in ISR 
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Abstract 

Using 13.3 fb^^ of e'^e~ collision data taken in the T(1S' — 45) region with the CLEO III detector 
at the CESR collider, a search has been made for the new resonance y(4260) recently reported by 
the BaBar Collaboration. The production of F(4260) in initial state radiation (ISR), and its decay 
into vr+TT" J/^, are confirmed. A good quality fit to our data is obtained with a single resonance. 
We determine M(y(4260)) = (4284l|^(stat) ± 4(syst)) MeV/c^, r(y(4260)) = (73+^^(stat) ± 
5(syst)) MeV/c2, and ree(i"(4260)) x ^(F (4260) ^ tt+tt^J/V') = (8.9l|?(stat) ± 1.8(syst)) eV/c^. 
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The BaBar Collaboration |l| has reported the observation of a broad structure in 
the TT^TT~J/il) mass spectrum around 4.26 GeV/c^ in initial state radiation (ISR) events, 
e"'"e~ — i> 7isRvr+7r~ J/-?/;. Data with a total luminosity of 233 fb~^ were taken at and 
near the peak of the T(4S') resonance. The enhancement, labeled y(4260), was fit- 
ted as a single resonance with a signal/background ratio of approximately 2/1, number 
of events N{Y) = 125 ± 13, mass M(F) = (4259 ± 8(stat) j:^(syst)) MeV/c^, width 
T(Y) = (88 ± 23(stat)l^(syst)) MeV/c^, and TeeiY) x B{Y n+n-J/ij) = (5.5 ± 
1.0(stat)lo;7(syst)) eV/c^. The observation of y(4260) in ISR leads to the conclusion 
that it is a vector with J^*^ = 1 , with additional evidence provided by the obser- 
vation that the dipion mass spectrum peaks at high mass, consistent with a S-wave 
phase space model. Recently, BaBar has reported that no evidence for y(4260) exci- 
tation is found in the reaction e^e~ 7isrPP) ^-nd obtained a ratio which is equivalent 
to a 90% confidence limit of Tee(Y) x B{Y pp) < 0.72 eV/c^ [2]. In the reaction 
J / ipTi^ -K" K~ they report a 3. la resonance signal and establish a 95% confidence 
limit of T{B- YK') x B{Y ^ tt+tt-J/?/') < 255 eV/c^ (sj. None of the theoretical 
calculations of charmonium spectra predict a vector state near 4260 MeV/ c^, and the mea- 
surements of the parameter R = (T(hadrons)/(T(/x"'"/i~) actually show a dip in this region, 
which is quite well fitted with the known charmonium vectors [4]. These facts make y(4260) 
an unexpected and provocative resonance, and many theoretical models have been proposed 
to explain its nature jl, 0] . This also makes it extremely important to make an independent 
experimental confirmation of the y(4260) resonance and its parameters. This paper reports 
on just such a confirmation. 

In a recent measurement of direct e~^e~ annihilation at ^/s = 4040, 4160, and 4260 
MeV/c^, CLEO R has confirmed the existence of an enhancement in the n~^7i J/ip cross 
section at 4260 MeV/c^ by measuring a cross section consistent with the peak cross section 
reported by BaBar. Similar enhancements were observed in n^ir^J/ijj and K^K~ J 
However, because the near-peak y(4260) data were limited to a single energy point, the 
mass and width of y(4260) could not be determined. In this paper, we present the results 
of the ISR production of y(4260) in 13.3 fb~^ of e^e~ data taken in the T(1S' — 45*) region 
with the CLEO detector at CESR, the e^e~ collider at Cornell University. Although our 
integrated luminosity is considerably less than that used by BaBar, the low background in 
our measurement enables us to make a statistically significant confirmation of y(4260), and 
an independent determination of its mass and width. 

The data used in the present investigations consist of 1.40, 1.77, 1.60, and 8.57 fb~^ taken 
at and in the vicinity of the T(IS'), T(2S'), T(3S'), and T(4S') resonances, respectively. The 
CLEO detector in its CLEO III |8| configuration was used to detect final state particles. 
The reaction studied was 

e+e~ 7iSR7r+7r"J/V^, J/^) e+e", 

The four charged particles in the final state, tt+tt" and were detected. It was required 
that both electrons, or at least one muon be identified in the barrel or end cap regions of the 
detector. The energy E deposited in the calorimeter, and the momentum p measured by the 
tracking chamber, were used to form E/p. It was required that for electrons 0.85 < E/p < 
1.1, and for muons E/p < 0.3. Electron candidates were required to have specific ionization 
energy loss, dE/dx, in the drift chamber consistent with the electron hypothesis within 3a. 
For muons a penetration depth of >3 hadronic interaction lengths in the muon chambers was 
required. Pions were identified using both the dE/dx information from the drift chamber 



3 



3850606-004 




-4 -2 2 4 6 8 101214 



MLs (GeV^/c^) 

FIG. 1: Distributions of data events (points) for M^^^g, as defined in the text. The peaks centered 
at M^jgg = and bounded by the dashed hnes correspond to the missing ISR photons. The shaded 
histograms show ip{2S) ISR signal Monte Carlo predictions normalized to the peak counts. 

and the information from the Ring Imaging Cherenkov (RICH) detector. The four charged 
particles were required to have net charge equal to zero. Further, very loose requirements 
were used, namely, total energy < 6.5 GeV, total transverse momentum, pt < 3.5 GeV/ c, 
and total energy of showers not matched with charged particles, -Eneut < 5 GeV. Lepton 
pairs were accepted in the range \M{l^l~) — M[JN))\ < 50 MeV/c^. The event selection 
criteria were optimized by using GEANT-based [9| F(4260) ISR signal Monte Carlo and 
J/ip side band yield as background, to maximize S'^/{S + B), where S and B refer to the 
signal and background event counts, respectively. The sideband was defined as 2.7 to 3.5 
GeV/c^, excluding the region 2.9 to 3.2 GeV/c^. 

The ISR events were identified by constructing the mass of the missing particle as 

M^.g, = (2Ebeam - E {n+ J / f - pin+ tt' J / i^f . 

The resulting M^j^^ distributions are shown in Fig. 1. The peaks centered at M^j^^ = 
correspond to the missing ISR photon and have little background under them. Above 
the peak in the spectrum for J/ip f^~^fJ'~, the background increases rapidly because of 
misidentified pions. We consider events with |M^jgg| < 1.7 GeV^/c'^ for J/ip e'^e~, and 
events with |M^igg| < 1.0 GeV^/ for J/ip ^ fJ'~^f^~, as ISR events. In order to improve the 
mass resolution, these events were fitted with the lepton tracks constrained to a common 
vertex, and the dilepton mass constrained to the J/ip mass, 3097 MeV/ c^. The distributions 
for the of the vertex and J/ip mass constraints, lepton and pion momenta, transverse 
momentum, and neutral energy for these selected events were found to be in good agreement 
with the corresponding distributions for the ip{2S) ISR signal Monte Carlo (more than 95% 
of the ISR events come from iIj{2S)). In Fig. 2 we show the M {71^71^ J ftp) distributions for 
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FIG. 2: Event distributions after kinematic fitting for the invariant mass M^tt^tt J/ip), with 
J/ip ^ e~^e~ and Evidence for the production of Y(4260) is present in both distributions. 



these events for J/ijj —>■ e^e and /i^/i . Both plots show clear enhancements around 4260 
MeV/c^ with very little background. This background is found to be in good agreement 
with that determined from the sidebands of J/ip. 

The mass resolution was determined by fitting the ISR signal Monte Carlo distributions 
for M {tt^ TT~ J / ip) for ip{2S) and y(4260), assuming zero intrinsic width for both. The in- 
strumental resolution, full-width at half-maximum, was thus determined to be 8.2 MeV/ 
for i){2S) and 10.6 MeV/c^ for F(4260). The data for ^(25) and Y(4260) enhancements 
were each fitted with a constant background plus the peak shape determined by the Monte 
Carlo simulation convoluted with the single resonance relativistic Breit-Wigner function 
which includes the phase space factor. The mass, width, peak area, and the background 
were kept free in the fits. The fit for the ip{2S) peak returned N{iIj{2S)) = 285 ± 17 counts 
and M{iIj(2S)) = (3685.70±0.24) MeV/ c\ which compares favorably with the known mass, 
M{ip{2S)) = (3686.093 ±0.034) MeV/c^ [ij. The fit for the F(4260) enhancement is illus- 
trated in Fig. 3. It results in M(F(4260)) = (4284^1^) MeV/c^ and r(y(4260)) = (73^1) 
MeV/c^. The number of counts in the fitted resonance is iV(y(4260)) = 13.6t3;9. The fit 
with a single Breit-Wigner has x^/d.o.f. = 1.00, a confidence level of 52%, and a signifi- 
cance of 5.4(7. The significance level for the resonance was obtained as a = y^— 2 ln(Lo / i^max) , 
where L^ax is the maximum likelihood value for the fit with resonance, and Lq is the likeli- 
hood value for the fit without resonance. 

The systematic errors in our measurements of M(y(4260)) and r(y(4260)) are much 
smaller than the statistical errors. An estimate of the systematic error in our mass mea- 
surement for y(4260) is provided by our result for M[il)(2S)). As mentioned earlier, our 
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FIG. 3: The M [tt'^ tt^ J / ip) distribution for the sum of J /tp — > e'^e~ and J /ip fi^/j,^ . The dotted 
hne denotes the fitted background, and the sohd curve denotes the fit with a single resonance. The 
semilog plot in the insert illustrates the tp{2S) excitation in the extended mass region. 
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FIG. 4: Dipion mass distribution for data events (points). The shaded histogram is the Monte 
Carlo prediction for M(7r''~7r~) based on a S'-wave phase space model. Neither the data nor the 
Monte Carlo distributions have been corrected for the experimental efficiency. 

measured mass M{'ip{2S)) differs from its known value by only (0.4 ± 0.2) MeV/c^. Ad- 
ditional contributions to the systematic uncertainty in M(y(4260)) and r(y(4260)) were 
estimated by fitting the y(4260) peak in the data by varying by it50% the resolution width 
returned by our ISR Monte Carlo sample. We have also studied the effect of variations in 
event selection cuts (M^^gg, pt, E^eut, for vertex and M{J/ip) constraints, lepton iden- 
tification criteria), and variations in peak shape and background parameterization. The 
maximum variation found was ±4 MeV/c^ in M(F(4260)) and ±5 MeV/c^ in r(r(4260)), 
which we therefore assign as systematic uncertainties. 
Thus our final results for the mass and width are 

M(F(4260)) = (4284+}^(stat) ±4(syst)) MeV/c^ 
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r(y(4260)) = (73l|(stat) ±5(syst)) MeV/c^ 
The quantity Tee{Y) x B{Y n^n^J/ip) is related to the observed counts N{Y) as 

N{Y)/B{J/^ ^ /+/-) 



ree(y) X B{Y ^ vr+vr- J/V^) 



C X e(F) X aisR(F) 



where £ is the integrated luminosity, e{Y) is the efficiency determined by the ISR signal 
Monte Carlo, and crisRiY) is the value of the ISR cross section per keV/c^ of di-electron 
width calculated using the theoretical formalism of Kuraev and Fadin 

From our sample of 4'{2S) ISR events, for which Teeii'i^S)) x B{ip{2S) n^n^J/tp) is 
known 13], we find that our ISR signal Monte Carlo efficiency e{ip{2S)) is overestimated by 



a factor of 1.5. We attribute this to the fact that because of the extreme forward nature of 
the ISR events the geometrical acceptance of our detector (| cos 6^1 <0.93) allows us to accept 
only approximately 15% of the phase space for il){2S), and we become extremely sensitive 
to the Monte Carlo modeling of the ISR production mechanism. For y(4260) events the 
geometrical acceptance of our detector is larger, approximately 30%, and we therefore expect 
that the Monte Carlo determined efficiency for the y(4260) events is more accurate than 
that for the ip{2S) events. We estimate the appropriate efficiency correction for the F (4260) 
events to be between 1.0 and 1.5, and we assume it to be 1.25±0.25. This is the dominant 
systematic uncertainty in our measurement of rf,f.{Y) x B{Y tt^tt' J/ip). Systematic 
uncertainties in luminosity, ctisr, and B{J/iIj — *• are much smaller. To further minimize 
the systematic uncertainty in our determination of ree(^) x B{Y —>■ -k^-k^ J/ip) we impose 
the additional requirement for its estimate that at least one lepton from the J /if) decay 
have |cos6'| <0.7. As a result, the data were refitted for the Y(4260) resonance. The fit 
returned N{Y) = S.ltH with a x^/d.o.f. = 0.98. With aisR(F) = 6.26 pb/keV, and 
e{Y) = 11.5%/(1.25 ± 0.25), our final result is 

ree(F) X B{Y ^ TT+TT- J/V^) 

= (8.9l|?(stat) ± 1.8(syst)) eV/c^. 

This is consistent with the BaBar result. 

Fig. 4 shows that the dipion mass of the y(4260) events is enhanced at high mass, as 
predicted by Monte Carlo generated distribution based on a phase space model modified by 
rescattering of pions which are in a relative S state [l2| . 

In summary, we have conffimed BaBar's observation of a new resonance, y(4260), 
produced in initial state radiation. We obtain resonance parameters in agreement with 
theirs. A good quality fit to our data is obtained with a single resonance. The observation 
of y(4260) production in ISR and the ^-wave phase space distribution of the dipion mass 
con&m that J^^(F(4260)) = 1~ [H, 0. The various theoretical models proposed for 
y(4260) have been discussed in detail in reviews by Eichten, Lane and Quigg [5], and 
Swanson [6|, but no clear understanding of the nature of y(4260) exists at present. A 
dedicated resonance scan and good measurements of various decays of F(4260) are needed 
to discriminate among the different theoretical models. 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
luminosity and running conditions. D. Cronin-Hennessy and A. Ryd thank the A. P. Sloan 
Foundation. This work was supported by the National Science Foundation, the U.S. De- 
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